Slow transport by continuous time quantum walks 
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Continuous time quantum walks (CTQW) do not necessarily perform better than their classical counterparts, 
the continuous time random walks (CTRW). For one special graph, where a recent analysis showed that in 
a particular direction of propagation the penetration of the graph is faster by CTQWs than by CTRWs, we 
demonstrate that in another direction of propagation the opposite is true; In this case a CTQW initially localized 
at one site displays a slow transport. We furthermore show that when the CTQWs initial condition is a totally 
symmetric superposition of states of equivalent sites, the transport gets to be much more rapid. 
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PACS numbers: 05.60.Gg,05.40.-a,03.67.-a 

The transfer of information over discrete structures (net- 
works) which are not necessarily regular lattices has become 
a topic of much interest in recent years. The problem is rel- 
evant to many distinct fields, such as polymer physics, solid 
state physics, biological physics and quantum computation, 
see Refs. |:jj,l2.l2l|4|,|j,|a] for reviews. In particular, quantum 
mechanics seems to allow a much faster transport than clas- 
sically possible. Thus, recent studies of quantum walks on 
graphs show that these often outperform their classical coun- 
terparts, i.e., in terms of the penetrability of the graph, for 
an overview see Ref. |4] and references therein. We recall 
that the extension of classical random walks to the quantum 
domain is not unique. There exist different variants of quan- 
tum walks, such as discrete |7] and continuous time ver- 
sions, which are not equivalent to each other. Here we focus 
on walks in continuous time. 

Walks occur over graphs which are collections of connected 
nodes. To each graph corresponds a discrete Laplace opera- 
tor (sometimes also called adjacency or connectivity matrix), 
A = (Aij). Here the non-diagonal elements Ay equal —1 if 
nodes i and j are connected by a bond and otherwise. The 
diagonal elements An equal the number of bonds which exit 
from node i, i.e., An equals the functionality fi of the node i. 

Classically, assuming the transmission rates of all bonds to 
be equal, say 7, the continuous-time random walk (CTRW) is 
governed by the master equation |6] 



dt 



Pj 



(*) =^2 T 3kPk(t), 



(1) 



where pj (t) is the probability to find at time t the walker at 
node j. T = (Tjk) is the transfer matrix of the walk, which 
is related to the adjacency matrix by T = — 7A. Equation Q 
is spatially discrete, but it also admits extensions to continu- 
ous spaces, e.g., leading to the disordered Lorentz gas model, 
which describes the dynamics of an electron through a disor- 
dered substrate |9]. 

We stick with the spatially discrete situation and let now the 
CTRW start from node k, i.e., we set Pfe(0) = Sjk- Denoting 
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FIG. 1 : (a) Graph consisting of two Cayley trees of generation G = 
3. (b) horizontal projection of the graph following Ref. ITch . (c) 
vertical projection of the same graph, (d) Vertical projection of a 
similar graph, obtained from two Cayley trees of general generation 
G, indicating the new nodes (clusters) and the dk, see text for details. 



by Pjk(t) the conditional probability of being at node j at time 
t when starting at node k at t = leads to 



(2) 



which is another way to write Eq.Q QlfJ]. Given the linearity 
of these equations, their solution involves a simple integration. 
For Eq.0 we have formally 



p jk (t) = {j\e"\k) 



(3) 



We now turn to one quantum mechanical extension of 
the problem, the so-called continuous-time quantum walk 
(CTQW). CTQWs are obtained by assuming the Hamiltonian 
of the system to be H = T iM flQIl . Then the basis vectors 
|fc) associated with the nodes k of the graph span the whole 
accessible Hilbert space. In this basis the Schrodinger equa- 
tion reads 
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i-|*) = H|fc), 



(4) 



2 



where we set H = 1. The transition amplitude aj k (t) from 
state \k) at time to state \j) at time t is then 

a jk (t) = (j\e- mt \k). (5) 

According to Eq.® the a.jk{t) obey 

d 



i—a 
dt 



(6) 



The inherent difference between Eq.Q and Eq.(|5j is, apart 
for the imaginary unit, the fact that classically . Pj k (t) = 1, 
whereas quantum mechanically J2j \ a jk(t)\ 2 = 1 holds. 

We turn now to the graph displayed in Fig.^a). The graph 
is obtained from two finite Cayley trees of generation G which 
have a common set of end nodes along the horizontal symme- 
try axis indicated in Fig. [Da), O3 Uflll . For the nodes on the 
axis as well as for the top and bottom nodes the connectivity 
is / = 2, whereas for all other nodes / = 3. 

The authors of Refs. d and CU have analyzed CTQWs 
over the graph given in Fig. 0a) > focussing on walks which 
start at the top node, and looking for the amplitude, Eq.(|5}, 
of being at the bottom node at time t. The problem can 
then be simplified by considering only states which are to- 
tally symmetric superpositions of states |fc) involving all the 
nodes fc in each row of Fig. 0a), as indicated schematically 
in Fig. 0b). The transport gets then mapped onto a one- 
dimensional CTQW m . 

Given that CTQWs obey time inversion, so that they never 
reach a limiting distribution, one uses the quantity 11111 



Xjk 



= lim — 



dt 



(*)| 2 



(7) 



to compare the efficiency of CTQWs to that of CTRWs. We 
will show in the following that the x?fe may depend strongly 
on the initial state. Now, as shown in 1 10], based on Eq.0, the 
CTQW's probability of being at the bottom node when start- 
ing at the top node is considerably larger than that of CTRWs. 

One legitimate question to ask now is: What happens if one 
considers on the same graph CTQWs which start at the left- 
most node and end at the rightmost node? As we proceed to 
show, it turns out that then the transport by CTQWs gets to be 
much slower than the transport by CTRWs. We start by fo- 
cussing on the full solution of Eq.(|6}, for which all the eigen- 
values and all the eigenvectors of T = H (or, equivalently, 
of A) are needed. For, say, the 22 nodes of Fig. 0a) we have 
to solve the eigenvalue problem for A (or T), which is a real 
and symmetric 22 x 22 matrix. This is a well-known problem, 
also of much interest in polymer physics 1311 . and many 
of the results obtained there can be used for our problem here. 

We recall first that the matrix A is non-negative definite. 
Then, for a structure like the one in Fig. 0a), A has exactly 
one vanishing eigenvalue, Ao = 0, the remaining eigenvalues 
being positive. Let A„ denote the nth eigenvalue of A and A 
the corresponding eigenvalue matrix. Furthermore, let Q de- 
note the matrix constructed from the orthonormalized eigen- 
vectors of A, so that A = QAQ -1 . Now the classical prob- 
ability is given by 

p jk (t) = (jlQe-^CT 1 ^). (8) 
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FIG. 2: Top: Graph obtained from two Cayley-trees of generation 
G = 2. Below: Probabilities for the CTQW, (a) and (b), and for 
the CTRW, (c) and (d), to be at node j after time t when starting 
at time from node 1 or from node 4. Left column, (a) and (c): 
Starting node is the top node 1. Right column, (b) and (d): Starting 
node is the leftmost node 4. The time is given in units of the inverse 
transmission rate 7" 1 , see text following Eq.0. 



For the quantum mechanical transition probability it follows 
that 



7r jk (t) = \a jk (t)\ 2 = \(j\Qe 



-it-yA 



Q 



(9) 



In order to determine numerically the corresponding eigen- 
values and eigenvectors of the matrix A for different graphs 
we have used the standard software package MAPLE 7. We 
start by considering the smaller graph, G = 2, given at the 
top of Fig. |2] The figures show the transition probabilities for 
CTQWs and CTRWs starting at the top node 1 (left column), 
which corresponds to the situation described in I1CB . or at the 
leftmost node 4 (right column). Remarkably, CTQWs start- 
ing at the top node reach the opposite node 10 very quickly, 
see Fig. |2j a), much quicker than expected from the CTRW be- 
havior, Fig.|3c). However, for walks starting at the leftmost 
node 4 and going to the rightmost node 7, the probabilities for 
the CTQWs, Fig. f£b), and the CTRWs, Fig. Hd), get to be 
comparable. Furthermore, the CTQWs' probabilities, 7^4 (i) 
and 7Ti i(f), of return to the starting node within the time in- 
terval depicted in Fig. [2] are much higher if the walks start at 
the leftmost node 4 of the graph instead of at the top node 1 . 
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FIG. 3: Ratios iTjk (i) /Pjk(t) for different directions of propagation, 
(a) top-bottom walk and (b) left-right walk over time t, see Fig|2|for 
units and details. 

On the other hand, for CTRWs there is not much difference 
between starting at the leftmost or at the top node, only that in 
the first case it just takes a little bit longer to reach an uniform 
distribution, compare Fig. |2jc) and Fig.|2jd). 

We now extend the time interval to t = 40 and compare the 
efficiency of the CTQW transport to the CTRW one. In Fig. [5] 
we plot for the top-bottom and the left-right walks the ratio of 
the quantum mechanical probabilities TTjk(t) to the classical 
ones Pjk(t). For top-bottom transport, depicted in Fig.^a), 
the plot turns out to be highly regular, reflecting the high sym- 
metry of the underlying graph in the vertical direction. For 
left-right transport the plot is less regular. Note the differ- 
ent scaling of the ordinates in the two parts of Fig. [5] which 
again stresses the preferential role played by the transport in 
top-bottom direction. 

In order to discuss what happens at even longer times, we 
proceed to evaluate for the CTQW the limiting distributions 
Xjk given by Eq. @. For CTRWs the limit is simple: All 
Pjk{t) tend to the same constant, which is the inverse of the to- 
tal number of nodes in the graph, no matter where the CTRWs 
start. For the CTQWs, however, this is not the case, as can al- 
ready be inferred from Figs. [2] and [3] Hence, we compute 
Xjk for the top-bottom and for the left-right cases separately. 
Therefore, for the G — 2 graph we compute the eigenvalues 
and eigenvectors of the respective 10 x 10 matrix A and the 
Xjk by using MAPLE. Having the appropriate eigenvalue ma- 
trix A and the matrix Q constructed from the orthonormalized 
eigenvectors we find with Eqs. @ and 10 and the LinearAl- 
gebra package of MAPLE that for the top-bottom case 

Xio.i = 0.2644 > i (10) 
whereas for the left-right case 

X7A = 0.0545 < ^. (11) 
Thus, the limiting CTRW-probability, 1/10, lies between X7A 




FIG. 4: Transition probabilities irjk{t) for a CTQW on the G = 3 
graph starting at (a) the leftmost node, TTjs(t), and (b) starting at the 
top node, ivji(t), for times t = 1, 5, 20, 80, and 160, see Fig|2|for 
units and details. 
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FIG. 5: Transition probability for a CTQW on the G = 2 graph start- 
ing at (a) the leftmost node, nj4,(t), and (b) at the top node, nji(t), 
for times t = 1, 5, 20, 80, and 160, see Fig|2|for units and details. 



and Xio.i- The top-bottom CTQW is more, the left-right 
CTQW less efficient than the corresponding CTRW. 

In order to better visualize that the top-bottom and left-right 
Xjk are very different, we show in Figs. |4] and [5] the quantum 
mechanical transition probabilities Ttjk (t) for all nodes j when 
starting (a) at the leftmost node 8 and (b) at the top node 1 ; in 
these figures the time is displayed parametrically. Figure|4]is 
for the G = 3 graph and Fig.|5]for the G = 2 graph. Now, 
even for the small graphs considered here, we find differences 
in the transition probabilities, which clearly depend on the ini- 
tial node. For the G = 3 graph consisting of 22 nodes, the 
CTQW starting at the top node 1 spreads out rapidly over the 
whole graph. After a very short time interval, there is a large 
probability to find the walker at the bottom node 22, see Fig. 
IHb). However, for the CTQW starting at the leftmost node 8, 
we have up to times t = 160 a high probability of finding it in 
the left half of the graph, see Fig.|4ja). Therefore, the prop- 
agation of the CTQW is strongly dependent on the starting 
node. For the smaller G — 2 graph of Fig. [5] which consists 
of 10 nodes, the effect is similar, but slightly less pronounced. 
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FIG. 6: Limiting probability Xjk for a CTQW on the (a) G = 2 
graph and on the (b) G = 3 graph. Starting points are the top node, 
k = 1 and the leftmost nodes / = 4 and ft = 8, respectively. 



We illustrate the situation at very long times in Fig.|6] where 
we display the limiting probabilities \jk for the G = 2 and 
the G = 3 graphs, see Eq.(0. For a CTQW starting at the top 
node 1 the limiting probability distribution has its maximum at 
the end nodes of the graphs, i.e. at nodes 1 and 10 for G = 2, 
and at nodes 1 and 22 for G = 3. For a CTQW starting at 
the leftmost node, k = 4 f or G = 2 and k = 8 f or G = 3, 
the limiting probability distribution shows a strong maximum 
around the starting node. 

Other initial conditions for the CTQW are, indeed, possible, 
especially when considering the high symmetry of the under- 
lying graphs. Note that, using for instance the site enumera- 
tion of Fig.0] a CTQW from node 8 to node 15 is equivalent 
to a CTQW from, say, node 10 to node 14. The graph's sym- 
metry suggests to collect groups of such nodes into clusters, 
while focussing on the transport from left to right. It is then 
natural to view the nodes 8, 9, 10, and 1 1 as belonging to the 
first cluster. The second cluster consists then of the nodes 4, 
5,16, and 17, all of which are directly connected by one bond 
to the nodes of the first cluster. The nodes 2 and 20 of the 
third cluster are all nodes directly connected by one bond to 
the nodes of the second cluster, while at the same time not 
belonging to the first cluster. In general, all the nodes of the 
(k + l)st cluster are connected by one bond to nodes of the 
fcth cluster and at the same time do not belong to the (k — l)st 
cluster. 

Let us denote the number of nodes in cluster k by dk- The 
transport occurs now from a cluster to the next, by which the 
original graph gets mapped onto a line in which one new node 
corresponds to a group of original nodes of the graph. For a 
new node at position k € [2, G] we find that dk = 2 G ~ fc+1 , 
the same being true for the mirror node value, i.e., dk — 
d 2 o+2-k- Note that for the end nodes d\ — d 2 c+i — 2 G_1 , 
the same holds for the nodes next to them. Moreover, for the 
middle node do+i = 2. 
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FIG. 7: (a) Transition probability TVji (t) for a CTQW between dif- 
ferent clusters j of the G — 3 graph. The CTQW starts at the first 
cluster, presented is the situation at times t = 1, 5, 10, 20, 40, 80, 
and 160, see Fig|2|for units and details, (b) Limiting probability Xji 
for a CTQW starting at the first cluster. 



We now focus on the transport via the states which are to- 
tally symmetric, normalized, linear state-combinations for all 
the original nodes in each cluster. Thus, for the fcth cluster, 
whose sites we denote by n, we have as a new state 



\a k ) 



1 



^dk , 

ft net 



(12) 



The CTQW is now determined by the new Hamiltonian 
H = 7 A, where the matrix elements of A are obtained from 
the new basis states \ak) and from the matrix A through 



Aj k = (oj-|A|afc) 



(13) 



Given the properties of A and the construction of the \ak), 
Ea. (ll2> . A is a real and symmetrical tridiagonal matrix, which 
implies a CTQW on a line. The diagonal elements of A are 
given by 

A kk = (a*|A|a fc ) = j- ^ («'|A|n) - /„ = f k , (14) 

where fk is the functionality of every node in the fcth cluster. 
For the sub- and super-diagonal elements of A we find 



ifc.fc+i 



Ak+i,k 
1 



y/dkdk- 



(a k \A\a k+ i) 



1 net 



\fd~kdk 



=05) 



fe+i 



where bf- is the number of bonds between the clusters k and 
k + 1. Now, except for the ends and the center of the graph, 
bk equals the maximum of the pair (dk, dfc+i). Between the 
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central node (c?g+i = 2) and its neighbors {do — c?g+2 = 2) 
the number of bonds is b G = &g+2 = 2. The number of 
bonds between the end node and its neighbor is b\ = &2G+1 = 
2di = 2 G . 

For the graph consisting of 22 original nodes the new ma- 
trix A is a tridiagonal 7x7 matrix, which can be readily di- 
agonalized. The advantage of the procedure is clear: the new 
matrix A depends on the number of clusters and grows with 
(2G+ 1), whereas the full adjacency matrix, A, grows with the 
total number of nodes in the graph, namely with (3 • 2 G — 2). 

From Ea.( ll2t the transition amplitude between the state 
|afc) at time and the state |<ij) at time t is given by 

a jk (t) = ( aj \e- tflt \a k ) = HQe^'Cr 1 !^) , (16) 

where A is the eigenvalue matrix and Q the matrix con- 
structed from the orthonormalized eigenvectors of the new 
matrix A. 



Now the quantum mechanical transition probabilities are 
given by 7fjfe (t) = |cr,fc(i)| 2 . FigureQa) shows the transition 
probabilities for CTQWs over clusters. Remarkably now, and 
similar to Fig.|3b), already in rather short periods of time such 
CTQWs move from one end cluster to the other. The limiting 
probability distribution, Xjk> which is depicted in Fig.Qb), 
also supports this finding. Note that Fig. 0b) again reflects 
the symmetry of the original graph. 

In conclusion, we have shown that CTQWs do not neces- 
sarily perform better than their CTRWs counterparts. By fo- 
cussing on a particular graph, we have shown that the pene- 
tration of such a graph by CTQWs can be better or worse than 
the one by CTRWs, depending on the initial state and on the 
propagation direction under scrutiny. 
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